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^ ■ Abstract 

\ The range of neutralino-proton cross sections for R-parity preserving super- 

gravity models with GUT scale unification of the gauge coupling constants 
is examined. The models considered here are mSUGRA, models with non 
universal soft breaking and D-brane models. It is found that the current dark 
matter detectors are sampling significant parts of the SUSY parameter space 
£T) '. and future detectors could sample almost the entire space. The special regions 

of parameter space that may be inaccessible to future detectors are seen to 
have a squark/gluino spectra beyond 1 TeV, but observable at the LHC. 
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While the existence of dark matter, which is now believed to make up about 30 % of the 
matter and energy of the universe, has been known for sometime, the nature of dark matter 
■ remains unresolved. Supersymmetric models with R parity invariance automatically posses a 
dark matter candidate, the lightest supersymmetric particle (LSP). We consider here models 
based on gravity mediated supergravity (SUGRA), where the LSP is generally the lightest 
neutralino (x?)- Neutralinos in the halo of the Milky Way might be directly detected by 
their scattering by terrestrial nuclear targets. Such scattering has a spin independent and a 
spin dependent part. For heavy nuclear targets the former dominates, and it is possible to 
extract then (to a good approximation) the neutralino -proton cross section, cr^o_ p . Current 
experiments (DAMA, CDMS, UKMDC) have sensitivity to halo Xi for 

a^_ p ~ 1 x l(T 6 pb (1) 

and future detectors (GENIUS, Cryoarray) plan to achieve a sensitivity of 

^o_ p ^(10- 9 -10- 10 )pb (2) 

We consider here two questions: (1) what part of the parameter space is being tested by 
current detectors, and (2) what is the smallest value of cr^o_ p the theory is predicting (i.e. 
how sensitive must detectors be to cover the full SUSY parameter space). The answer to 
these questions depends in part on the SUGRA model one is considering and also on what 
range of theoretical and input parameter one assumes. In the following, we examine three 
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models that have been considered in the literature based on grand unification of the gauge 
coupling constants at M G = 2 x 10 16 GeV: (1) Minimal super gravity GUT (mSUGRA) with 
universal soft breaking at Mq 0; (2) Nonuniversal soft breaking models with Higgs and third 
generation scalar masses at Mq p|, and D-brane models (based on type IIB orientifolds ||) 
which allow for nonuniversal gaugino masses and nonuniversal scalar masses at Mq @]. 

While each of the above models contain a number of unknown parameters, theories of 
this type can still make relevant predictions for two reasons: (i) they allow for radiative 
breaking of SU(2) x U(l) at the electroweak scale (giving a natural explanation of the Higgs 
mechanism), and (ii) along with calculating cr^o_ p , the theory can calculate the relic density 

of Xii i- e ^x? = Px°/Pc where p^o is the relic mass density of Xi an d p c = 3Hq/8ttGn {H 
is the Hubble constant and Gat is the Newton constant). Both of these greatly restrict the 
parameter space. In general one has fl^oh 2 ~ (/ i dx(a ann v }) _1 (where <T ann is the neutralino 
annihilation cross section in the early universe, v is the relative velocity, Xf — kTf/m^o, Tf 
is the freeze out temperature, (...) means thermal average and h = F /100 km s^Mpc" 1 ). 
The fact that these conditions can be satisfied for reasonable parts of the SUSY parameter 
space represents a significant success of the SUGRA models. 

In the following we will assume Hq = (70±10)km s _1 Mpc _1 and matter (m) and baryonic 
(b) relic densities of Vt m = 0.3 ± 0.1 and tt b = 0.05. Thus Vt^oh 2 = 0.12 ± 0.05. The 
calculations given below allow for a 2a spread, i.e. we take || 

0.02 < Q^oh 2 < 0.25. (3) 

It is clear that when the MAP and Planck sattelites determine the cosmological parameters 
accurately, the SUGRA dark matter predictions will be greatly sharpened. 

I. CALCULATIONAL DETAILS 

In order to get reasonably accurate results, it is necessary to include a number of theo- 
retical corrections in the analysis. We list here the main ones used in the calculations below: 
(i) In relating the theory at Mq to phenomena at the electroweak scale, the two loop gauge 
and one loop Yukawa renormalization group equations (RGE) are used, iterating to get a 
consistent SUSY spectrum, (ii) QCD RGE corrections are further included below the SUSY 
breaking scale for contributions involving light quarks, (iii) A careful analysis of the light 
Higgs mass rrih is necessary (including two loop and pole mass corrections) as the current 
LEP limits impact sensitively on the relic density analysis for tan/3 < 5. (iv) L-R mixing 
terms are included in the sfermion (mass) 2 matrices since they produce important effects 
for large tan/3 in the third generation, (v) One loop corrections are included to m& and m T 
which are again important for large tan/3, (vi) The experimental bounds on the b — > sj 
decay put significant constraints on the SUSY parameter space and theoretical calculations 
here include the leading order (LO) and approximate NLO corrections. We have not in the 
following imposed b — t (or t — b — r) Yukawa unification or proton decay constraints as these 
depend sensitively on unknown post-GUT physics. For example, such constraints do not 
naturally occur in the string models where SU(5) (or S'O(IO)) gauge symmetry is broken by 
Wilson lines at M G (even though grand unification of the gauge coupling constants at Mq 
for such string models is still required). 
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All the above corrections are under theoretical control except for the b — > s'y analysis 
where a full NLO calculations has not been done. (We expect that while the full analysis 
might modify the regions of parameter space excluded by the b — > 57 experimental con- 
straint, the minimum and maximum values of cr^o_ p would probably not be significantly 
changed.) The analysis of cr^o_ p , taking into account the above theoretical corrections has 
now been carried out by several groups obtaing results in general agreement |6HT2j. These 
results are presented below. 

Accelerator bounds significantly limit the SUSY parameter space. In the following we 
assume the LEP bounds pi m h > 104(100) GeV for tan/3=3(5) and m + > 104(100) GeV. 



(For tan/3 > 5, the rrih bounds do not produce a significant constraint fll3fl .) For b — > 57 we 



assume an allowed range of 2a from the CLEO data [14 



1.8 x KT 4 < B(B X sl ) < 4.5 x 10~ 4 



(4) 



The Tevatron gives a bound of nig > 270 GeV( for niq = nig) [jT5] . 

Theory allows one to calculate the x^-quark cross section and we follow the analysis of 
|TS| to convert this to Xi ~ P scattering. For this one needs the ir — N a term, 



VttN = -K + m d )(p\uu + dd\p), 



(5) 



Co = o'ttN — (jn-u + rrid) (p\ss\p) and the quark mass ratio r = m s j \l/2){m u + nid). We use 
here a n N = 65MeV, from recent analyses |T7| , [T8"|j based on new n — N scattering data, cr = 30 
MeV p| and r= 24.4 ± 1.5 @. Older n — N data gave a wN 

p by a factor of about 3. 



45 GeV [211, which if used 



would reduce the overall cr^o 



II. MSUGRA MODEL 

We consider first the mSUGRA model where the most complete analysis has been done. 
mSUGRA depends on four parameters and one sign: mo (universal scalar mass at Mq), mi/2 
(universal gaugino mass at Mq), Aq (universal cubic soft breaking mass, tan/3 = (H 2 ) / (Hi) 
(where (H 2 ,i) gives rise to (up, down) quark masses) and where \i is the Higgs mixing 

parameter in the superpotential, = [iHiH 2 ). One conventionally restricts the range of 
these parameters by "naturalness" conditions and in the following we assume m < 1 TeV, 
mi/2 < 600 GeV (corresponding to m g < 1.5 TeV, m^o < 240 GeV), [Ao/m | < 5, and 
2< tan(3 < 50. Large tan/3 is of interest since SO (10) models imply tan/? > 40 and also 
(T^o_ p increases with tan/?. u^o- p decreases with m 1 / 2 for large m^, and thus if one were 
to increase the bound on m 1 / 2 to 1 TeV {rrig < 2.5 TeV), the cross section would drop by a 
factor of 2-3. 

The maximum cr^o_ p arise then for large tan/? and small m\/ 2 . This can be seen in Fig.l 
where (a^o_ p ) mauX is plotted vs. m^o for tan/3=20, 30, 40 and 50. Fig. 2 shows ^o/i 2 for 
tan/3 = 30 when the cross section takes on its maximum value. Current detectors obeying 
Eq (1) are then sampling the parameter space for large tan /3, small m^o and small ^^oh 2 
i.e 

tan(3 ~ 25, m^o ~ 90GeV, Vt^h 2 ~ 0.1 (6) 
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To discuss the minimum cross section, it is convenient to consider first m^o ~ 150 GeV 

Al 

( m i/2 < 350) where no coannihilation occurs. The minimum cross section occurs for small 
tan/3. From Fig. 3 one sees 

a ~o_ p £ 4 x 10~ 9 pb; m^> ~ 140GeV (7) 

which would be accessible to detectors that are currently being planned (e.g. GENIUS). 

For larger m^o, i.e. ~ 150 the phenomena of coannihilation can occur in the relic 

density analysis since the light stau, f 1; (and also e#, /t^) can become degenerate with the 
Xi- The relic density constraint can then be satisfied in narrow corridor of mo of width 
Am ~ 25 GeV, the value of m increasing as m 1 / 2 increases ||. Since m and mi/ 2 increase 
as one progresses up the corridor, cr^o_ p will generally decrease. 

We consider first the case fi > One finds in general that cr^o_ p also decreases as 



A increases. Fig. 4 shows c^o_ p in the domain of large A and for two values of tan/3. One 
sees that the smaller tan/3 still gives the lower cross section, though the difference is mostly 
neutralized at larger mm. (For large tan/3, mo also becomes large to satisfy the relic density 
constraint i.e m = 700 GeV for tan/9=40, mi/ 2 = 600 GeV.) We have in general for this 
regime 

a^o_ p ~ 1 x 10~ 9 pb; formi/ 2 < 600GeV,/i > 0, A < 4m 1/2 . (8) 

This is still within the sensitivity range of proposed detectors. 

When fi is negative an "accidental" cancellation can occur in part of the parameter 
space in the coannihilation region which can greatly reduce 0£o_ p ||. This can be seen in 
Fig. 5, where starting with small tan/? the cross section decreases, leading to a minimum 
at about tan/3=10, and then increases again for larger tan/3. At the minimum one has 
(T^o_ p = 1 x 10~ 12 for tan/3=10 and myz = 600 GeV. More generally one has 

a k o_ p < 1 x 10~ 10 pb (9) 

for the parameter domain when 4 ~ tan/3 ~ 20, mi/2 ~ 450GeV(mg ~ l.lTeV), fi < 0. 
In this domain, cr^o_ p would not be accessible to any of the currently planned detectors. 
However, mSUGRA also then predicts that this could happen only when the gluino and 
squarks have masses greater than 1 TeV (and for only a restricted region of tan/3) a result 
that could be verified at the LHC. 



III. NONUNIVERSAL SUGRA MODELS 

In the discussion of SUGRA models with nonuniversal soft breaking, universality for the 
first two generations of squark and slepton masses at Mq is usually maintained to suppress 
flavor changing neutral currents. One allows, however, the Higgs and third generation squark 
and slepton masses to become nonuniversal. The masses may then be parametrized at Mq 
as follows: 

m£=mg(l + <5i); = m 2 {l + 5 2 ); (10) 

m q L = m l( l + h)] rn t 2 R = mjj(l + 5 4 ); rn T 2 R = + <J B ); 
m 4 =m o( 1 + J 6); m^ = mg(l + 5 7 ). 
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where qi = (£l, II = (^t, Tl)- Here m is the universal mass of the first two generations 
and 5{ are the Higgs and third generation deviations from universality. In the following we 
will restrict the 5i to obey 

- 1 < Si < 1. (11) 

The lower limit on Si is necessary to prevent tachyonic behavior, and the condition Si < 1 is 
taken as a reasonable upper bound. We also maintain gauge and gaugino mass unification 
at Mq. 

While these models contain a large number of new parameters, their effects on <J^o_ p can 
be charcterized approximately by the signs of 8%. ..84 flU] , Thus the choice of S3, 64, S\ < 



and 82 > can greatly increase c^Q-p; and the reverse choice can reduce cr^o_ p ( though by 
a much lesser amount). The possible increase is shown in Fig.6 for tan/? = 7, fi > where 
nonuniversal soft breaking increases cr^o_ p by a factor of 10-100 compared to the universal 
case. Thus it is possible for detectors to probe regions of smaller tan/3 with nonuniversal 
breaking, and detectors obeying Eq. (1) can probe part of the parameter space for tan/? as 
low as tanft ~ 4. 

The minimum cross section occurs (as in mSUGRA) at the lowest tan/3 and at the largest 
mi/2 i.e. in the coannihilation region. We limit ourselves here to the case where only the 
Higgs masses are nonuniversal i.e. 7^ 0. One finds then results similar to mSUGRA i.e. 
cr^o_ p ~ 10~ 9 pb for jj, > 0, mi/2 < 600 GeV. For fi < 0, there can again be a cancellation of 
matrix elements reducing the cross section to 10~ 12 pb when TO1/2 = 600 GeV in a restricted 
part of the parameter space when tan/3 ~10 ||23|| . 



IV. D-BRANE MODELS 

Recent advances in string theory have stimulated again the construction of string inspired 
phenomenologically viable models. One type, based on type IIB orientifolds puts the 
Standard Model on sets of 5 branes. These models are of interest in that they can lead to a 
pattern of soft breaking different from what can naturally arise in conventional supergravity 
GUTs. An interesting example |4j] has the following soft breaking pattern at Mq: 

fhi = m 3 = — A = \/3cos8b9ie~ iai m3/2 (12) 
m 2 = VScos9 b (l - ©i) 1/2 m 3/2 

where and are the gaugino masses, and 

m i2 = (! - 3/2sin 2 9 b )ml /2 forq L , l L , H h H 2 (13) 
m\ = (1 - 3sin 2 9 b )m 2 3/2 foTU R , d R , e R . 

Thus the SU (2) doublets are all degenerate at M G but different from the singlets. We note 
Eq. (13) implies 9 b <0.615. 

We consider first the case with no CP violating phases in the soft breaking sector, i.e. 
«i=0 and the fi parameter real. In general cr^o_ p is a decreasing function of 9 b since the 
squark and slepton masses increase as 9 b decreases. This is illustrated in Fig. 7 where cr^o_ p 
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is plotted for tan/3=10, m 3 / 2 = 175 GeV with 9b = 0.5 (upper curve) and 9b = 0.2 (lower 
curve). One has that current detectors obeying Eq. (1) are sensitive to the model in the 
domain when tan/3 ~ 15. However m^o is quite small when tan/? is close to the minimum 
value. 

The lower bound on cr^o_ p should occur when 9b and tan/3 are small. Fig. 8 shows the 

minimum cross section for m^o < 300 GeV, fi > 0. One sees <J^o_ p ~ 1CT 9 pb. As in 
mSUGRA, a cancellation of matrix elements can occur for < 0, allowing a smaller cross 
section in the domain 7 ~tan/3 ~ 30 and reaching a minimum of ~ 10~ 12 pb for m 3 / 2 = 600 

GeV at tan/3 ~ 12. These low cross sections again imply very heavy gluinos i.e m g ~ 950 
GeV. The D-brane model also possesses an interesting new region of coannihilation for 
©! ~ 0.8 when the light chargino Xi can become degenerate with which has not been 
analyzed here. 

CP violating phases produce contributions to the electric dipole moments of the electron 
d e , and neutron d n . The bound on d e is d e < 4.3 x 1CT 27 ecm at 95 % C.L. [p^j , and one 
must constrain the parameter space to satisfy this. This generally reduces c"^o_ p by a factor 



of 2-3 ||12|| . However models with CP violating phases require a great deal of fine tuning 
unless tan/3 < 3-5, so one would expect in this case cr^o_ p to be quite small. 



V. SUMMARY 

We have examined here the neutralino-proton cross section for a number of SUGRA type 
models. In all the models considered, there are regions of parameter space with Xi~P cross 
sections of the size that could be observed with current detectors. Thus with the sensitivity 
of Eq. (1), detectors would be sampling regions of the parameter space for mSUGRA where 
tan/3 ~ 25, m^o ~ 90GeV and fi^o/i 2 ~ ^-l- Nonuniversal models can have larger cross 

sections and so detectors could sample down to tan/3 ~ 4, while for the D-brane models 
considered, detectors could sample down to tan/3 ~ 15. 

The minimum cross sections these models predict are considerably below current sensi- 
tivity. Thus for mSUGRA one finds for /i > 

a^o_ p ~ 1 x 10" 9 pb form 1/2 < 600GeV, /i > 0. (14) 

where mi/ 2 = 600 GeV corresponds to m g = 1.5 TeV, m^o = 240 TeV. This is still in the 
range that would be accessible to detectors being planned (such as GENIUS or Cryoarray). 
For jj, < 0, a cancellation can occur in certain regions of parameter space allowing the cross 
sections to fall below Eq. (14), Thus 

a^_ p < 1 x l(T 10 pb for 4 ~ tan/3 ~ 20, [i < 0, m 1/2 ~ 450GeV (15) 

and reaching a minimum of cr^o_ p = 1 x 1CT 12 pb for tan/3=10, m!/ 2 =600 GeV, /i < 0. This 
domain would appear not to be accessible to future planned detectors. Since mi/ 2 =450 GeV 
corresponds to m g = 1.1 TeV, this region of parameter space would imply a gluino squark 
spectrum at the LHC above 1 TeV. 

The above results holds for the mSUGRA model. While a full analysis of coannihilation 
has not been carried out for the nonuniversal and D-brane models, results similar to the above 
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hold for these over large regions of parameter space. Thus for nonuniversal Higgs masses 
and for the D-brane model one finds cr^-p ~ 10~ 9 pb for /i > 0, while a cancellelation allows 
o"^o_ p to fall to 10~ 12 pb for \i < at tan/3 ~ 10 (with again a gluino/squark mass spectrum 
in the TeV domain). 

In the above we have limited m\/2 to m\/2 < 600 GeV (m§ < 1.5 TeV). Increasing mi/2 
lowers the cross section. Thus allowing mi/ 2 = 1 TeV {rrig = 2.5 TeV) would reduce the 
minimum cross section by a factor of about 2-3. Also using the earlier value of a n N = 45 
MeV 



21 rather than the most recent value ex. 



ttN 



65 MeV 



n],|l8[ would reduce cr^,o_ p by a 



factor of 3. 
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FIGURES 
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FIG. 1. (a ■yO 15 J max 

for mSUGRA obtained by varying Aq over the parameter space for 
tan/3 = 20, 30, 40, and 50[11]. The relic density constraint, Eq.(3) has been imposed. 
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FIG. 2. ^Xi^ 2 f° r m SUGRA when (a^o_ p ) takes on its maximum value for tan/? = 30 [11]. 
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FIG. 3. (a ^o_p) m i n for Tri-^o < 140 GeV for mSUGRA obtained by varying Aq for tan/3 = 3 
with the relic density constraint, Eq.(3)[25]. 
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FIG. 4. (<T^o_ p ) for mSUGRA in the coannihilation region for tan/3 = 40 (upper curve) and 
tan/3 = 3 (lower curve), Aq = 4m 1 / 2) A 4 > 0[25]. 
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FIG. 5. (cr%o_ p ) for mSUGRA and fi < for (from top to bottom on right) tan/3=20, 5 and 
10. Note that for tan/3 > 10, the curves terminate at the left due to the b — > constraint [25]. 
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FIG. 6. Universal mSUGRA (lower curve) and nonuniversal soft breaking (upper curve) for 
tan/3 = 7, fj, > 0. The nonuniversal curve has <5~3, 62, 5i < 82 > 0[11]. 
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FIG. 7. a x°~P f° r ^ e D-brane model for tan/3 = 10, m 3 / 2 = 175 GeV, \i > 0. Upper curve is 
9 b =0.5, lower curve b =O.2[12]. 
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